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ABSTRACT 


We study the production of a Higgs boson recoiling from a massless invisible system in e + e - 
collisions. This is a quite distinctive signature that can arise when the Higgs boson is produced 
in association with a massless dark photon, which can happen in BSM scenarios foreseeing an 
extra unbroken U( 1) gauge group. Dark photons can indeed acquire effective couplings to the 
Higgs boson as occurs in models recently proposed to generate exponentially-spread Yukawa 
couplings. We analyze the signal and corresponding backgrounds for H —> bb, and estimate 
ILC and FCC-ee sensitivities in a model-independent way. 



1 Introduction 


The LHC discovery of the Higgs-boson resonance at 125 GeV [T] has definitely strengthened our 
confidence in the Higgs mechanism as the origin of the electroweak (EW) symmetry breaking 
(EWSB) and fermion mass generation [2]. All present data are well consistent with the Standard 
Model (SM) expectations for the Higgs boson [3]-|4|, although there is still room for potential 
New Physics (NP) effects. 

At the same time, the absence of any NP signal at the LHC Run I is causing considerable 
concern about the applicability of the naturalness criteria. The latter would require new phe¬ 
nomena at the TeV scale to stabilize the SM Higgs scalar potential against potentially large 
radiative corrections coming from NP energy thresholds. The SM is also facing the Flavor 
problem, which is related to the unexplained huge hierarchy in the fermion mass spectrum or, 
analogously, in the Higgs Yukawa couplings. 

On another front, there is increasing evidence from astrophysical and cosmological observa¬ 
tions of the existence of Dark Matter (DM) in the Universe p, ;6], which is not predicted in the 
SM. Independently of naturalness criteria, there might then be NP above the EW scale that 
explains DM and, in some cases, could be tested at the LHC. 

A common origin for DM and Flavor is also conceivable. One can postulate the existence 
of a hidden (dark) sector, where all these issues are addressed, which is composed of new fields 
that are SM singlets. The Higgs boson can then act as a portal to the dark sector J7]. The 
Flavor and EWSB structures are indeed restricted to the Higgs couplings and mass, and are 
not related to other SM couplings. On the other hand, NP could well affect the Higgs-boson 
characteristics by smaller amounts than the present LHC sensitivity in Higgs-boson data. 

In this connection, in jBJ a new paradigm has been proposed to generate exponentially- 
spread Yukawa couplings from gauge quantum numbers in the dark sector. In this class of 
models, the Flavor and chiral symmetry breaking (ChSB) take place in a dark sector, and 
are transmitted by radiative corrections to the observable sector through Higgs-portal type of 
interactions. The Yukawa couplings arise radiatively as effective couplings at low energy. The 
hidden sector consists of (stable) massive dark fermions (that are SM singlets and potential 
DM candidates), and a massless dark photon, the gauge boson of an unbroken U(1)f gauge 
group in the dark sector. Chiral symmetry is spontaneously broken, and dark fermions obtain 
non-vanishing (flavor-dependent) masses via a non-perturbative mechanism involving U(1)f 
gauge interactions. The resulting chiral-symmetry and flavor breaking in the dark sector is 
then transferred to the Yukawa-coupling sector at one-loop via scalar messenger fields that 
are charged under both SM and U(1)f gauge interactions. A similar framework has also been 
explored in jHl HUJ, although no unbroken U( 1) gauge sector is introduced in that case. 

The new unbroken U (1) gauge group (and the corresponding massless dark photon) is a cru¬ 
cial dynamical component of the model in [8], but is also a common feature of various theories of 
new physics, including models with gauge-symmetry breaking of compact gauge groups, string- 
theory motivated phenomenological models, and models of interacting dark matter mm- 
It is indeed conceivable that a hidden sector contains an extra long-range force. Remarkably, 
being massless, an on-shell dark photon can be fully decoupled from the SM quark and lepton 
sector at any order in perturbation theory m- This is not true for a massive dark-photon, 
due to a potential tree-level mixing with the photon field. Most of present astrophysical and 
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Figure 1: Higgs decays -H —)• 77 , Z 7 , 77 , via mediator loops. 


accelerator constraints [13] apply to massive dark-photon couplings, and can be evaded in a 
massless dark-photon scenario. This allows for potentially large dark-photon couplings to the 
dark sector, that might also lead to observable new signatures at colliders [ 8 j. 

The Higgs boson can interact with dark photons radiatively. In the framework proposed in 
[ 8 ], this occurs at one loop by the exchange of scalar messenger fields (Figure [l]). 

As a consequence, the Higgs boson can act as a portal toward the dark sector, giving rise 
to new Higgs-boson decays such as pH] 

H 77 , Z 7 , 77 , (1) 

where the symbols 7 and 7 stand for the usual QED photon and dark photon, respectively, and 
Z is the neutral vector boson. The corresponding decay rates can in principle be large, even for 
very heavy messenger fields. As in the H —> 7 7 , Z 7 , gg decays in the SM, the non-decoupling 
Higgs properties guarantee non-vanishing decay widths even in the large-mass limit for particles 
exchanged in the loop, provided the virtual (messenger) fields carry the same SU{2)l quantum 
numbers of quarks and leptons. 

Being fully decoupled at tree level from the SM sector, a single massless dark photon will give 
rise in the Higgs final state to same amount of missing energy and missing momentum, while 
the two dark-photon channel will contribute to the invisible Higgs rate. Extra contributions to 
the widths of the SM channels H 77 , Z'y^gg are also expected in general. 

The H —> 7 7 decay gives rise to a new spectacular signature at the LHC in 7 + fix final 
states, with a photon plus missing transverse energy fi T resonating at the Higgs mass. In pH] , 
a parton-level study shows that the LHC Run-1 data set could be sensitive to BR (H —>- 77 ) 
values as low as 0.5%, while a minimal-model prediction for BR (H —> 77 ) can be as large 
as 5%. 

The aim of the present study is to analyze the phenomenological implications of the Higgs 
effective couplings to dark photons at future e + e _ colliders usmsun]. Apart from the new 
signatures corresponding to the Higgs-boson exotic decays in Eq. (|TJ) (that we do not address 
here), effective H 77 and H'yZ interactions involving dark photons will give rise to final states 
with a Higgs boson and a dark photon, 

e + e _ —> H 7 , (2) 

via s-channel exchange of either a photon or a Z vector boson (Figure [2]). This channel, 
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Figure 2: Diagrams for e + e —> F/ 7 . 


although kinematically similar to the SM one-loop channel e + e _ — > Hy |18j . gives rise to a 
completely new signature, since the final massless 7 goes undetected. 

We will focus on the bb 7 final state corresponding to the main Higgs decay channel H —y bb, 
although even more rare Higgs decays will be of relevance in the clean e + e“ environment mm- 
The e + e~ —y Hy final state will then be characterized by an unbalanced bb system resonating 
at the Higgs mass rrin, the dark photon 7 giving rise to “monochromatic” missing energy $ 
and momentum (for fixed initial c.m. collision energy a/s). Contrary to what occurs in 
the main irreducible SM bbvv background, at parton level the invariant mass of the invisible 
system M miss = {$" — p 2 ) 1 ^ 2 vanishes. This feature will provide a crucial handle for background 
suppression. 

Since the messenger fields are expected to be quite heavy with respect to the characteristic 
energy of the e + e“ — y H'y process, the H 77 and HZ 7 vertices can be considered as effective 
interactions, and described by a model-independent parametrization na. The ratio of the HyZ 
and Hyy couplings will in general depend on the spin and the SM gauge-group representation 
of the new particles running in the loop. For simplicity, we will focus here on scenarios where 
the H'yZ vertex is induced by scalar messenger fields in the SU(2) L x SU( 3) c fundamental 
representation [Hj, which gives a definite prediction for the H'yZ and H 77 coupling ratio. 

The paper is organized as follows. In Section 2, we provide a model-independent parametriza¬ 
tion of the effective couplings controlling the Higgs exotic decays H —> 77 , Z 7 , 77 , and the 
SM-like decays H —s- 77 , Zy, and express the relevant Higgs BR’s in terms of the model- 
independent coefficients. In Section 3, we study the sensitivity of future e + e _ colliders to the 
e + e _ —>• H 7 associated production by analyzing the signal and corresponding backgrounds. In 
Section 4, we discuss the NP model in [5] that aims to solve the Flavor hierarchy problem. We 
also present the corresponding predictions for the Higgs-dark-photon effective couplings, and 
for the Higgs branching ratios (BR’s) relative to the decays H —> 77 , and H —> Zy. Finally, 
our conclusions are discussed in Section 5. In the Appendix, we describe some U(l)p coupling 
properties of the model in [ 8 j, that are needed to discuss its phenomenological consequences. 

2 Effective dark-photon couplings to the Higgs boson 

We now introduce the dark-photon effective couplings to the Higgs boson that enter the e + e~ —>■ 
Hy cross section. I 11 general, Higgs-dark-photon effective couplings can arise at one loop due 
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Figure 3: Effective coupling approximation for the vertices H 77 , HZ 7 . 


to the exchange of messenger fields that are charged under both the SM and the U(1)f gauge 
groups (Figure |3]). In case the messenger masses are much larger than both run and y/s, one 
can use the effective theory approximation. The corresponding effective Lagrangian £^ sgs can 
be split as 

£rfr &gb = £dp h + £sm h > (3) 

where £dp h contains the dark-photon effective interactions with the Higgs boson, while £ SMh 
presents the extra (that is messenger-induced) contributions to the SM Higgs effective interac¬ 
tions with two photons, one photon and a Z , and two gluons. 

By retaining only the relevant low-energy operators, £dp h can be expressed in terms of 
dimensionless (real) coefficients Cij (with i , j = 7 , 7 , Z,g) as 

£ = (4) 

71 \ V V V J 

where a is the SM fine structure constant, and 7 ^, Z^ U) 7 ^ are the held strentghs of photon, 
Z boson, and dark photon, respectively ( 7 ^ = d^A v — d u A^ for the photon held A^). Then, 
£sm h can be written as 

£sm h = + ^Z^H) + — C A^LG ai ' lJ G^H, (5) 

where as is the SM strong coupling constant, G a ^ u stands for the gluon held strength, and a 
sum over the color index a is understood. 

As usual, the Cij coefficients in Eqs.(|4])-([5j can be computed in the complete theory by 
evaluating one-loop amplitudes for relevant physical processes, and by matching them with the 
corresponding results obtained at tree level via the effective Lagrangian in Eq. (J3J) . In partic¬ 
ular, in order to express the coefficients C 7 ^, C 77 in Eq.Q in terms of the fundamental 

parameters of the model, one can match the tree-level widths, based on the parametrization in 
Eq.Q, for the Higgs decays H 77 , H —>■ Z 7 , H —» 77 , respectively, with the corresponding 
one-loop results computed in the full model (as sketched in Figure [3]). This will be discussed 
in Section 4, after introducing a particular NP framework. 

On the other hand, one can perform a phenomenological study of the e + e _ —> Hj process 
just on the basis of the model-independent parametrization in Eq.Q, which we will do in the 
next Section. 
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Figure 4 : BR’s for H 77 and H —y Z7, normalized to the SM values for BR (H —» 77) and 
BR(iJ —» Z7), respectively, versus C 77 and Cz 7 - 


Before proceeding, we connect the basic C l3 coefficients in Eq.Q to the corresponding 
H —» * j decay widths. The H —>■ 77 width has been computed in ra. and, taking into account 
the parametrization in Eq. Q , one has 


r(H ^ 77) = 


m H a 


! |C 77 , 


87 T 3 U 2 


T(H gg ) = 


<4|CWI 2 (^ 2 ~ 1) 

47 T 3 U 2 


( 6 ) 


where N c = 3 and T(H —> gg) is understood to be inclusive in gluons final states. Analogous 
results can be obtained for the H —> 77, H —> Z 7, H —» 77 widths replacing |C' 77 | 2 by 2 |C' 77 | 2 , 
|C'^ 7 | 2 , 21C 77 1 2 respectively. 

Figure [ 4 ] shows the branching ratios for H —>■ 77 and H —> Z7, normalized to the SM 
BR(// —» 77) and BR(i/ —> Z7), respectively, versus the corresponding Cij coefficients. The 
Cij ranges shown in the plot include values well allowed by the model described in Section 4 . 
One can then get for the Higgs decays into a dark photon an enhancement factor ( 9 ( 10 ) with 
respect to the SM Higgs decays where the dark photon is replaced by a photon. This makes 
the corresponding phenomenology quite relevant for both LHC and future-collider studies. 

It is also useful to express the BR’s for H —$■ 77, 77, 77 as a function of the relative exotic 
contribution to the H —>■ ij decay width, as the ratio 


r ij = 


rm 

1 ij 

pSM 

77 


( 7 ) 


with T™- generically indicating the pure messenger contribution to H —> ij , with i,j = 7,7. 
Then, one obtain the following model-independent parametrization of the H —y 77, 77, 77 


6 








BR’s as functions of r, 7 m 



= BR™ 


tdtjSM 

— nn , 77 

y 

— -fc>/l 77 


1 + ryyERtof ’ 


1 + r V/ BR™ ’ 

(1 + Xy/r^Y 

1 + ryyBR* M ’ 


( 8 ) 


where x = ±1 parametrizes the relative sign of the SM and exotic amplitudes, and BR^ stands 
forBR 

Analogously, the relative deviation for the H —s- gg decay width will be defined as 


r 
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(9) 


3 Sensitivity study for e + e — H 7 

We focus now on the 7 production in association with a Higgs boson in e + e~ collisions. The 
e + e” —>• Hy total cross section versus a/s is shown in figure [ 5 ] for three different coupling 
assumptions: C 77 = 1 ,Czj = 0 (blue line); C 77 = 0, Czy = 1 (green line); C 77 = 1, CU 7 = 
0.79 (7 77 (red line). The coupling ratio Czy/ C 77 ~ 0.79 is typical for scenarios where the HyZ 
vertex is induced by scalar messenger fields in the SU(2) L x SU( 3) c fundamental representation 
(see Section 4). The corresponding cross sections at y/s ~ 1 TeV (relevant for linear colliders 
at larger collision energy) are 43 ab, 15 ab, 55 ab, respectively. Cross sections can be easily 
extrapolated to coupling set-up obtained just by globally rescaling these set of couplings. 

The e + e _ —> Hy cross sections grow with c.m. energy thanks to the nature of the dimension- 
five operators in the effective Lagrangian in Eq. Q. Hence, at constant integrated luminosity, 
higher-energy colliders will have a higher potential, since the dominant background is expected 
to scale down with energy as 1/s. On the other hand, lower y/s may allow larger integrated 
luminosity, as is the case of the e + e _ Future Circular Collider (FCC-ee) (also called TLEP) [T7j . 
where an integrated luminosity of 10 ab -1 is expected at a/s = 240 GeV. At linear colliders, 
either 1LC [15] or CLIC [16], one typically foresees integrated luminosities of a few hundreds 
fb _1 in the initial energy of yfs ~ 250 GeV or 350 GeV, and a few ab -1 at the larger —yfs 
stages [21]. Here, we assume the minimal energy setup of y/s = 240 GeV that is relevant for 
Higgs-boson studies, and study the sensitivity to e + e~ —> H'y production versus integrated 
luminosities foreseen at different machines. 

Using the effective Lagrangian Eq. Q implemented by FeynRules [22], we have gener¬ 
ated e + e~ —» H'y —> bby events with MadGraph5_aMC@NLO [23], and passed these events to 
PYTHIA to account for parton showering, and hadronization. We checked that the inclusion 
of effects from initial state radiation, that tends to degrade the c.m. energy in a circular 
e + e _ colliders, would moderately affect the results of the present analysis. We neither include 
beamstrahlung effects that can be of some relevance at linear colliders. We account for finite 
detector resolution by applying the jet-energy smearing a(E)/E = 30 Vo/VE, which is typical 
for ILC-kind of detectors [24] . 
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Figure 5: Total e + e~ —y if 7 cross section as a function of the c.m. collision energy, for different 
sets of effective couplings. 


The dark photon escapes the experimental apparatus undetected, and the final signal con¬ 
sists of two 6 quarks and large missing energy and momentum ft. In our simulation we 
reconstruct the missing momentum from the vector sum of all visible final-state particle mo¬ 
menta, after applying PYTHIA. In a lepton collider a H —» bb final state is not swamped by 
large QCD backgrounds, as occurs in hadronic collisions. Therefore, bb final states are the 
best channel to search for i /7 production, thanks to the H —)• bb large rate. After showering 
and hadronization, we reconstruct jets (and 6 -jets) according to the basic PYTHIA jet-cone 
algorithm, assuming a quite large cone aperture Rj = 1.5, which optimizes mass reconstruction 
[25]. The basic event selection is given by 

p b T > 20 GeV , \rjb\ < 2.5 , A R(bb) > 0.4 , $ > 40 GeV, (10) 

where A R(bb) = A if + A </> 2 is the angular distance between two 6 -tagged jets. We assume a 
6 -tagging efficiency of 80%, and a corresponding fake 6 -jet rejection factor of 100 for light jets. 

The main SM background for the bb + $ final state is given by the uubb production. This 
includes the on-shell processes ZZ —>■ uubb, ZH —> uubh, which give an almost monochromatic 
66 -pair system (similarly to the signal), and the vector boson fusion channel Hvv. A subdom¬ 
inant contribution comes from vvqq (mostly from on-shell Z pairs), where both light jets are 
mis-tagged as 6 jets. 

There are two kinematical variables that turn out to be particularly efficient in separating 
the signal from the background. First, we introduce the variable Mjj as the invariant mass of 
the two jets with largest px- This is directly connected to the 6 -pair invariant mass, and can 
be used to pinpoint events with 6 -quarks coming from Higgs decays, out of the smaller-Af^- 
events arising from Z —> bb. There is anyway part of the vvbb background that goes through 
the ZH production resonating at Mjj ~ ran, just as in the signal case. This is well illustrated 
by Figure [ 6 j where the normalized invariant-mass distributions of the 66 system are compared 
for signal and backgrounds. Second, we introduce the missing-mass variable Ai miss , defined as 
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Figure 6: Invariant-mass Mjj distributions for the two jets with largest Pt for the the sig¬ 
nal (solid line) and the two backgrounds vvbb (dashed line) and vvqq (dot-dashed line) after 
PYTHIA showering, hadronization, and jet-energy resolution effect. All distributions are nor¬ 
malized to 1. 


M miss = (11) 

where $ = sfs — Jj) ^visible and 0 = — P visible are the final-state missing energy and missing 
three-momentum vector, respectively (the sum over visible objects here includes both jets 
and lower-energy particles escaping jet reconstruction). The M miss variable is expected to 
approximately vanish in the partonic description of e + e _ —» Hj, corresponding to the massless 
invisible dark photon. A cut on M miss then proves to be remarkably efficient in further separating 
the signal from the main background, where M miss mostly matches an invisible zf-boson decaying 
into neutrinos. 

The Mmi SR spectrum of the signal and background processes are compared in Figure [7j 
after applying PYTHIA showering, jet reconstruction and jet-energy resolution effects on top 
of parton-lcvel simulation (right panel). The parton-level spectrum, shown in the left panel 
of the same figure, shows a distinct peak at M miss ~ 0 for the signal, and at M miss ~ M z for 
the background processes. No energy-resolution effect has been applied in the latter case, and 
the smearing of the peaks is just due to the presence of neutrinos from b decays, and to the 
possible off-shellness of the vv system in the background. Applying the parton showering, jet 
reconstruction and energy-resolution effects (as in the right panel of the figure) degrades the 
Mmis s spectrum of the signal quite a lot, shifting the peak away from zero and smearing it. 
Hence, an optimal detector resolution would be particularly crucial in this analysis. 

On the basis of the Mjj and M miss distributions in Figures [6] and [7J we set a suitable event 
selection. We require the invariant mass Mjj to be within 10% of the Mjj peak value of the 
simulated signal events, and then impose the missing mass to be below 40 GeV. The latter 
cuts make the vvqq background negligible. The vvbb background can still be slightly reduced 


9 


















M miss (GeV) M miss (GeV) 

Figure 7: M miss distributions for the the signal (solid line) and the backgrounds uvbb (dashed 
line) and vvqq (dot-dashed line) after parton level simulation (left) and after PYTHIA show¬ 
ering, hadronization and jet energy resolution effect (right). All distributions are normalized 
to 1. 


Process 

Cross section (fb) 

Acceptance after cuts (%) 

i /7 (C Z y = 0 ) 

10.1 x lO " 3 C'A 

17.3 

H'y (C/y = 0) 

4.8 x 10 ' 3 C\^ 

17.3 

ZZ 7 (C /7 = 0-79 C/y) 

13.8 x 10 ~ 3 CY 

17.3 

SM vvbb 

115. 

0.08 


Table 1: Cross sections (in fb) and corresponding acceptances after kinematical cuts on 

signal and SM background at a/s =240 GeV. Applied cuts include the initial event selection in 
Eq. (10), Mjj to be within 10% of the Mjj peak value of signal events, M miss < 40 GeV, and 
$ < 100 GeV. Cross sections include BR (H —* bb) ~ 0.58. 


after these cuts by making a further cut on the missing energy $. The $ spectrum is shown in 
Figure [8] for the signal and background events satisfying the previous Mjj and M miss cuts. Both 
the signal and background distributions peak at around the same value, with the background 
moderately shifted to larger $ values. Thus we require the missing energy to be below 100 
GeV. Including the initial event selection criteria, we altogether impose that the missing energy 
satisfies the condition 40 GeV < $ < 100 GeV. 

Table [l] shows the cross sections and the acceptances for the signal and the vvbb background 
after applying the cut-flow just described, for a fs =240 GeV. The signal acceptance is practi¬ 
cally insensitive to to a change in the relative contribution of the C and Cz 7 couplings. The 
corresponding acceptance for the vvqq background is negligible. 

On the basis of the Table [T| acceptances, we can work out the expected sensitivity to the 
signal for given values of the C 77 , Cz^ couplings. As usual, we define the signal significance as 
S/y/S + B, being S and B the event numbers for signal and background, respectively. Figure [ 9 ] 
shows the integrated luminosity needed to make a 5 a observation of the i /7 production in 
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Figure 8: $ distributions for the the signal (solid line) and the background vvbb (dashed line) 
after PYTHIA showering, hadronization and jet-energy resolution effect, and after applying the 
Mu and cuts described in the text. All distributions are normalized to 1. 



Figure 9: The 5 cr-sensitivity bounds for e + e^ —y H'y as a function of the integrated luminosity 
at a/s = 240 GeV, when C = 0 (green), Cz*/ = 0 (blue), and Cz y = 0.79 [with C shown 
on the vertical axis] (red). 
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Figure 10: Signal significance S/\/S + B for e + e~ —)■ H 7 as a function of the couplings 
Cjj,Cz^, for (7 77 = 0 (green), C z ^ = 0 (blue), and C z ~ f = 0.79 C 77 [with C 77 shown on the 
horizontal axis] (red), for an integrated luminosity of 10 ab -1 , at y/s =240 GeV. The horizontal 
gray lines show the 5 cr-discovery bound, and the 2cr (~ 95% C.L. exclusion) level. 


e + e^ collisions at y/s =240 GeV, for any given value of the C 77 , C z ~/ couplings (shown on the 
y- axis) when C 77 = 0 (green line), C z ^ = 0 (blue line) and C z ^ = 0.79C 77 (red line). 

For an integrated luminosity of 10 ab -1 at y/s = 240 GeV (a typical value for FCC-ee), 
Figure 10 shows the signal significance as a function of the couplings, with the same color 
convention as in Figure [9| The horizontal gray lines show the 5 cr-discovery bound on couplings, 
and the 2 cr level approximating the 95% confidence-level exclusion. 

Then, at 95% C.L., one can exclude the ranges C 77 > 1.9 (for C z ^ = 0), C z ^ > 2.7 (for 
C 77 = 0), and C 77 > 1.6 (for C z ^ = 0.79 C 77 ). The interval C 77 > 1.9 corresponds to a Higgs 
BR into 77 that is more than 3 times the SM BR (H —> 77 ), while C ^ 7 >2.7 corresponds to a 
Higgs BR into Z 7 that is more than 9 times the SM BR (H —>■ Z 7 ). 

The corresponding sensitivities on the C 77 and Cz 7 couplings at the ILC (foreseeing an initial 
y/s =250 GeV phase, with a typical integrated luminosity of 250 fb -1 ) can be estimated from 
Figure |9j and are about a factor 3 lower than the FCC-ee ones. The latter match sensitivities 
on the corresponding BR (H —> 77 ) and BR (H —> Z 7 ) that are smaller than the FCC-ee ones 
by an order of magnitude. 


4 A model of Flavor with Dark Photons 

In this section, we review the main aspects of the model proposed in jS], that provides a theoreti¬ 
cal framework for the effective description given by the Lagrangian in Eq. ([ 3 ]) . Correspondingly, 
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we will obtain predictions for the Higgs BR’s for H —y 77 and H —> Z 7 (and relevant effective 
couplings) in terms of the fundamental parameters of the model. 


4.1 The Lagragian 

The basic assumptions of the model introduced in [ 8 j are the following. For each SM fermion 
there exists a dark fermion with same flavor in a hidden sector, that is an heavier fermion replica 
which is a singlet under the SM gauge group. The chiral- and flavor-symmetry spontaneus 
breaking is realized in the dark sector as described in the following, and communicated to 
the SM fermions via renormalizable and flavor universal interactions mediated by messenger 
fields. The Yukawa couplings are radiatively generated, preserving approximately the same 
flavor hierarchy structure of the dark-fermion sector. A similar approach has been proposed in 
[91 HO], although in this case the dynamics responsible of the dark-fermion spectrum has not 
been discussed. 

A simple choice for the messenger sector consists of a set of scalar messenger fields that, 
due to gauge invariance, have to be charged under the SM gauge group with the same quantum 
numbers as quarks and leptons. The relevant Lagrangian is | 8 ] 


C 


— fY=o 4. r 
^SM ' *-"mes 


+ £dSi 


( 12 ) 


where C\^ stands for the SM Lagrangian without tree-level Higgs Yukawa couplings, C mes is 
the Lagrangian containing the messenger sector with its couplings to the SM and dark fields, and 
Cos is the dark-sector Lagrangian including the dynamics responsible of the flavor hierarchy. 

Dark fermions acquire an exponentially-spread mass spectrum from non-perturbative dy¬ 
namical effects. In fact, the Lagrangian Cos, containing dark fermions and a dark photon, is 
given by [ 8 ], 


Cos = iJ2(Q UiV ^Q Ui + Q DiV ^Q Di ) 

i 

- \f^ + (13) 

where Q Ul , Q Di are the dark fermion fields, partners in the hidden sector of the up (U t ) and 
down (Di) quarks, D M = d^ + igQA^ is the covariant derivative related to the U(1)f gauge field 
which is associated to the dark-photon, with F^a the corresponding field strength, and Q 
the charge operator acting on the dark fermion fields. Cos can be extended to include also the 
SM leptonic sector in a straightforward way. 

The last term in Cos , involving only the U(l)p gauge sector, corresponds to the so-called 
Lee-Wick term. It is a higher-derivative term, and the A scale can be interpreted as the mass 
of the associated massive ghostQ This term can trigger ChSB, and generate a mass spectrum 
non-perturbatively [26], by means of the Nambu-Jona-Lasinio mechanism [27]. The following 

According to the Lee-Wick argument J28J-[25], the presence of a massive ghost field in the spectrum does 
not spoil unitarity, provided the massive ghost has a finite decay width, which is automatically satisfied in the 
present scenario. 
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Dirac-fermion mass spectrum can be induced on the true vacuum m 


Mn. = A exp 


2vr 


1 

3b(A )qf + 4 


( 14 ) 


where g* is the U(l)p charge of the Dirac fermion, and d(A) the corresponding fine structure 
constant at the A energy scale. This solution is manifestely non-perturbative as can be seen 
from the a dependence in Mq.. 

The Lagrangian £ mes in Eq. (12) contains the messenger scalar fields, 


Cr 


= c + c 

*^mes 1 mes i 


(15) 


where £^ es is the kinetic Lagrangian for the messenger fields interacting with the SM gauge 
bosons, while £ l mes contains the messenger interactions with the dark fermions and the Higgs 
boson, which give rise to the effective Yukawa couplings. 

The SM quark quantum numbers set the minimal matter in the messenger sector, which is 
given by 

• 2 Nf complex scalar SU(2)l doublets: S L R and S '/ 1 , 

• 2 Nf complex scalar SU(2) L singlets: S R and S R , 


• one real SU{2)l x U( 1 )y singlet scalar: So, 

/ gUi(Di) \ 

where S R i(r>i) = ^ = 3, and i (= 1, 2, 3) stands for the flavor index for three 

fermion generations. The and S^ ,L>i fields carry the SM quark quantum numbers, and 

the labels L, R corresponds to the chirality of the associated SM fermions. They couple to the 
EW gauge bosons and to the gluons, as do squarks in the minimal supersymmetric extensions 
of the SM. Note that a minimal flavor violation would require this Lagrangian to be invariant 
under SU(Np), where Np = 2 Nf is the number of flavors. 

The messenger mass structure can be described by four free universal mass terms in both 
the S^' R and S'jj Ti sectors. Note that an even more minimal hypothesis of a common scalar 
mass for the L and R scalar sectors is also phenomenologically acceptable. 

The Lagrangian CJ mes for the messenger interactions with quarks and SM Higgs boson is 


£ 


i 

mes 


N f 


N f 


= 9l 


9r 


E fees] Y + E fees] Y + 


i =1 
N f 


i =1 
N f 


E [fees] ss + E [fee?] sg + 


2=1 


2=1 


+ \ s S„ (i/'SSfe*' + ff'SSfe*’) + h.c. + ViSij). 


(16) 


where S 0 is a real singlet scalar, V ( S 0 ) its potential, q l L , u l R , d 1 r , stand for the SM fermions, and 
H is the SM Higgs doublet, with H = Contractions on color indices are understood. 
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Fields 

Spin 

SU{ 2) l 

U{l)r 

St/(3) c 

V( 1 )F 

S°l 

0 

1/2 

1/3 

3 

-Qm 

S U L 

0 

1/2 

1/3 

3 

-Qui 

Sr* 

0 

0 

-2/3 

3 

-Qm 

S°n 

0 

0 

4/3 

3 

-Qui 

Q D > 

1/2 

0 

0 

0 

(hx 

Q Ui 

1/2 

0 

0 

0 

Qui 

So 

0 

0 

0 

0 

0 


Tabic 2: Spin and gauge quantum numbers for the strongly-interacting messenger fields and 
for dark fermions. U(1)f is the gauge symmetry in the dark sector. 


The two constants gi and g R are flavor-universal free parameters that are allowed to be in the 
perturbative region g LR < 1. 

The Lagrangian for the interaction of messenger scalars with the SM gauge bosons follows 
from the universal properties of gauge interactions. We stress that the messenger fields carry 
the same U(1)f charges as the corresponding dark fermions. 

The messenger and dark-fermion quantum numbers are shown in Table [2] (detailed in 0 )- 
Finally, after the spontaneous symmetry breaking of the discrete H —y —H parity symmetry 
that prevents tree-level Yukawa couplings, all Yukawa couplings are generated at one loop, and 
are finite at any order in perturbation theory mum)- Assuming almost degenerate diagonal 
messenger masses in the L and R sectors, one finds, from the one-loop vertex computation, 
that the effective Yukawa coupling associated to the quark i is [8] 


Y i = Y 0 (xi) exp 


2tt 


3« (A )q} 


2 1 ’ 


(17) 


where the dark-fermion masses Mq. have been replaced by Eq.(14), and the one-loop function 
Y 0 (xi) is given by 


= (fit) (ty) <*(*>• 


(18) 


with fig = A s{S ), and Xi = Mqjm 2 . Also, m is the average mass of the messenger fields 
running in the loop, and Cq{x) — (1 — x (1 — logic))/(l — x) 2 ^ 


4.2 BR (H -A 77 ) predictions 


The C 77 , C 77 , and C 77 coefficients entering the effective Lagrangian in Eq.([3| have been 


com¬ 


puted in [Hj, as a function of the basic parameters of the model described in Section 4.1 

In the m 2 L ~ fh 2 R approximation of degenerate messenger masses in the left and right-handed 
sectors, corresponding to the mixing angle 9 = n/4 (see [8] for notations), the flavor-universal 


2 The above results hold for diagonal Yukawa couplings. They can be easily generalized to include the 
Cabibbo-Kobayashi-Maskawa (CKM) matrix entering the charged weak interactions as explained in [9]. 
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messenger mass matrix can be expressed in terms of two parameters, the average messenger 
mass m 2 = ( m 2 L + m 2 R )/2, and the mixing parameter £ = A 2 jm 2 . We then define the universal 
mixing parameters £ g and £, corresponding to the messenger mixing parameters in the quark 
and lepton sectors, respectively. Note that, in the effective theory approximation, the Higgs 
and Z masses can be both set to zero in loop functions, when terms 0(in 2 H /m 2 ) are negligible. 

Then, one findtj^] 


C - 


Cyy 


iAJt 

12 1 -£| 


a 


E 

i=q,l 


2 ’ 


Q' z —' 12 1 — £ 

i=q,l 


2 ’ 
i 


C z y = 


c 77 = cl 




R\ C 2 


Z 1 


i=q,l 


12 1 - <£1 


2 ’ 


SM 

77 






i=q,l 


3FU-0 ’ 


Cz 7 = Cl 


SM 

77 


i + £«: 




i=q,l 


c gg = 


r ySM I i 4- 

99 1 3F q (i-e q ) r 


Zl 


£ 2 


3FC-G) 


where C™ 


|F, Cgg 1 = RF q , and the constants Rq \ 2 are given by 


(19) 


F* = 31V c ( e 2 +e 2 ), 

3 

R\ = N c 'y ^ (c-uQui T a D q D .^ , 
2=1 
3 


R 


q 

2 


N c + q n,) > 

i= 1 


/? — 3 p 2 
n 0 — o , 


-^1 — & e ^ ) , 

2=1 

R 2 = J2 « + q l) > 


2—1 


( 20 ) 


with = 2/3, e D = —1/3, and e s = —1, the electric charges for up-, down-quarks, and charged 
leptons, respectively. F and F q are the usual SM loop factor given by 

F = F w (/ 3 w ) + F f , F F = Y J N c Q)F f {t3 f ), F q = Y, F f(Pf), ( 21 ) 

/ / 

with N c = 1(3) for leptons (quarks) respectively, flw = 4m^/m#, (3f = 4 m 2 /m R , and 

F w (x) = 2 + 3x + 3x(2 -x)f(x), F f (x) — — 2x (1 + (1 — x)f(x)) , ( 22 ) 

3 Note that, due to the Bose statistics of the messenger fields, the relative sign with respect to the SM 

contribution in the H -4 gg amplitude is predicted to be negative. Analogously, there is a negative relative sign 

with respect to the SM fermion contribution to the H —> 77 and H -4 7 Z amplitudes. 
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l-\-y/l—X 
l — y/l—X 


, forx < 1. Including 


where f(x) = arcsin 2 ^], for x > 1, and f{x) = — | 
only the W and top-quark loops in F, we get \F\ ~ 6.5, \F q \ ~ 1.25 for m h = 125 GeV. We will 
elaborate on the constants in Eq. (19) (parametrizing the ratio of the messenger couplings 
to the Z and 7 ) at the end of this sectior^J 

One can see that in the Higgs couplings in Eq.(19), there is a clear non-decoupling effect, 
since the Cij coefficients do not vanish when (m 2 , A) —> 00 , provided the ratio A/m 2 is finite. 


The ratios r, 


J. LLK, ± I IJ (* , J — I , I ) 1 UV.1111V.U >11 11 

tion in Eq. (J 8 |, and r gg , defined in Eq. (|9]), are then given by 


7 , 7 ), defined in Eq. ([7]), entering the model-independent BR’s parametriza- 


r 77 — 


r 77 



r 77 — 


r 99 = 



X q F 2 

~F~ 


(23) 


(24) 


where the extra factor 2 in r 77 comes from statistics and 


X l(q) = 


e 


1(9) 


3F(1 




(25) 


i(g)' 


with R.q'j 2 defined in Eqs. (20) 


The strength of the exotic contribution to FI —> 77 is directly controlled by two mixing 
parameters, £ q and On the other hand, the H —>■ gg depends only on £ q , and can be 

constrained at the LHC by measuring Higgs production rates. 

It is useful to connect the messenger-loop impact on the Hgg vertex expressed by r gg [as 
defined in Eq. (J9j)] with the usual k g anomalous coupling of the Hgg interaction which enters 

/Rgg. Present 


the relation C gg = k g C g ^. By Eq. (19), it is straightforward to see that k g ~ 1 


data constraints k g at 68 % of C.L. to be in the ranges k g = l.OOlg^g (ATLAS Collaboration 
[3]), and k g = 0.76/q [3 (CMS Collaboration [3]). In the following, we assume r gg to be in the 
range 0 < r gg < 0.4. 


In Figure 11, we show BR(if —>■ 77 ) as a function of the U(1)f fine structure constant 
evaluated at the average messenger mass, a(rh), for a few values of the ratios r gg and r 77 . We 
assume the charge normalization qu 3 = 1 and = 1. The dashed blue lines correspond to the 
condition 


^BR sm 77 < BR 77 < 2BR 4 S 77 


> SM 


(26) 


with BR™ = 2.28 x 10 3 . The red dots correspond to fixed BR 77 ~ BR inv values, where 


BR inv is the Higgs invisible-decay BR. Note that when colored messengers contribute to the 
L /77 effective coupling, BR 77 depends also on the sign of the U(l)p charges, which are free 


parameters. Correspondingly, in Figure 11, we show cases in which the U(l)p charge sign in 


4 We will neglect the a running from the m scale to the characteristic low-energy scale entering the dark- 

photon vertex in H —> 77,77, Zg. 
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Figure 11: Predictions for BR (H — y 77 ) as functions of the U(1)f fine structure constant a, for 
different values of BR inv , r 77 and r gg in the full model; r 77 satisfies the constraint r 77 > 1.9 r gg 
(see text). 


the quark sector gives either destructive or constructive interferences with the electromagnetic 
charges. The maximum value of a ~ 0.18 corresponds to d(A) ~ 1 (see Eq.(34) in the Appendix 
for details). 

In the upper-left plot in Figure 11, we show the RR 77 predictions for r gg = 0. The allowed 
RR 77 values are at most about 1%, and one has £>R 77 > BR^ only for r 77 > 0.1, corresponding 
to a quite large mixing parameter (£; > 0.82) in the leptonic messenger sector. The upper-right 
plot shows the case r gg ~ 0.1, assuming constructive U(l)p charge interferences, corresponding 
to sign (qujq Di ) = sign (q Vi /q Ei ) = —1. In this plot, the minimum r 77 value ~ 0.2 comes from 
the constraint r 77 > 1.9 r gg , arising from Eqs. (24)-(25) as a result of the color and EW quantum 
numbers of the different amplitudes. RR 77 ahove 2% are allowed in this case, corresponding 
to the range a ~ 0.10 — 0.18, with r 77 ~ 0.4. 

The lower-right plot matches the largest allowed contribution from colored messengers ( r gg ~ 
0.4, with r 77 > 0.76), which corresponds to a mixing parameter £ q ~ 0.88, while the lower-left 
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corresponds to r gg ~ 0.2 with r 77 > 0.38. The left and right plots correspond to the destructive 
and constructive effects of the U(1)f charges, that is sign {quJq Di ) = sign (q Ui /q Ei ) = 1 , and 
sign (qirJq Di ) = sign {q v Jq Ei ) = —1, respectively. For large colored messenger contributions, 
the actual BR 77 value dramatically depends on the charge signs. In the case of constructive 
interferences, BR 77 reaches the 3%—4% level, while in the destructive case it is always below 
0.3%. 

Summing up, large BR 77 values up to 1%—4% are possible in this scenario. They correspond 
to the a range expected from naturalness arguments applied to U(l)p charges (see the Appendix 
for more details). This BR 77 range is equivalent to C 77 values up to about 5 in the Lagrangian 
in Eq. Q, as shown in Figure [4j As we have seen, in Section 3, these C 77 values could be 
well inside the domain of sensitivity of the e + e _ —> H 7 searches at future e + e _ colliders (cf. 
Figures [9] - [To|. 

Before closing this section, we elaborate on the general connection between the HZ 7 and 
i /77 couplings, induced at one loop. This will depend on the spin and the SM gauge-group 
representation of the particles running in the loop. On the other hand, here the ffZy and 
if 77 vertices are induced by scalar messenger fields in the SU(2)p x SU(3) C fundamental 
representation (cf. Table [2]), which gives a definite prediction for the HZ 7 and f/Vy coupling 
ratio, and for the f?j^ terms entering the C '^ 7 and C ^ 7 effective couplings in Eq.(19). 

Then we can now motivate the i ?^ 7 ~ 0.79 scenario that we considered in the ""eve - —>• if 7 
study of Section 3. Rz 1 can be defined (in absence of scalar mixing between the doublet and 
singlet scalars circulating in loops) by the relation 


Czx — Rz 7 C y x , 

with A" = 7 , 7 , for the Higgs effective couplings in Eqs. 0 - 0 - 
The Z-boson coupling to a scalar particle i is given by 



(!-§)- sin 2 °w 
sin 9 W cos 9\y 


(27) 


(28) 


times the photon coupling to the same scalar, where Y l and Q l are the hypercharge and the 
electric charge of the scalar, and sin 9w is the Weinberg angle. 

If the scalars in the loop are SM particle partners (as happens in SUSY or in the model 
in 0 ), they will share the quantum numbers of the left- and right-handed SM fermions. 

Then, for right-handed fermion partners (Y R = Q R ), and for left-handed electron partners 
( Y eL = — | and Q eL = —1), one has, respectively, 


R 


R _ 

z 7 


Mi 9 W 
cos 9 W 


~ -0.55, 


JD e L _ 2 

sin 9 


sin 2 9 W 
W COS 9\y 


~ 0.64. 


(29) 


The average Z-to -7 
leptonic scalars is then 


coupling ratio for a pair of mass-degenerate 



pi? 1 p«i 
^Z'y + ^Zj 

2 


~ 0.045. 


right- and left-handed 


(30) 


19 







For left-handed up- and down-type squarks ( Y UL = \Q UL and Y dL = —\Q dL , respectively), one 


has instead 


- | - sin 2 6 

jdu l _ _ 

Zl sin 6 W cos 6 


w 


~ 1 . 2 , 


sin 2 9 


w 


R dh = - 

Zl sin 0\y cos 6\y 


~ 3.0 


(31) 


'W '■'Ub U W 

The average contribution from a mass-degenerate pair of right- and left-handed up and down 
squarks is then R Zl = (R Zl + R^)/2 ~ 0.34, and R d Zl = (R Zl + R z L ry )/2 ~ 1.23, respectively. 
Assuming that also the up- and down-type scalars are mass degenerate, the net result from a 
squark doublet is then 


Dlt I JDd 

_ R Z 1 + U Zl 
K Zl - 9 


~ 0.79. 


(32) 


The same pattern for the R Zl constants can be obtained in the model in [8], in the approxi¬ 
mation of degenerate colored messenger scalars. 

In Section 3, we include the case C Z ^ = R Z/y C = 0.79 Cyy among the benchmarks for the 
analysis of the e + e _ —y H'y potential, corresponding to negligible leptonic contributions in the 
messenger loops. 


5 Conclusions 

Hidden sectors extending the SM theory can include an extra unbroken U( 1) gauge symmetry. 
The corresponding gauge boson, a massless dark vector boson, can couple to the Higgs boson 
through renormalizable interactions involving scalar messengers, giving rise to effective H'yj, 
HZ 7 , and H 77 couplings. Since a massless dark photon is not revealed by collider detectors, 
the latter can be probed at the LHC and future colliders via the search for exotic Higgs decays 
into a photon or a Z boson plus missing transverse energy, and the determination of the invisible 
Higgs decay width. 

Another way to investigate the possible existence of the Hjj and HZ') couplings is the 
production at future e + e^ colliders of a Higgs boson associated to a dark photon. The corre¬ 
sponding signature is very distinctive, since there is no irreducible SM background where the 
Higgs decay products are accompanied by a massless invisible system. 

After introducing an effective Lagrangian description of the new Higgs interactions, we 
studied the potential of the e + e~ —y H 7 production for probing the corresponding T /77 and 
HZ 7 couplings, C and C Z ~ ( . A most useful strategy for enhancing the S/B ratio turns out 
to be a selection on the small values of the invisible-system invariant mass. We found that, 
at a fs = 240 GeV, with the integrated luminosity foreseen at the FCC-ee (~ 10 ab _1 ), one 
can exclude at 95% C.L. the ranges CW > 1.9 (for C z ~ f = 0), C z ^ > 2.7 (for = 0), and 
Cyy > 1.6 (for C Z y = 0.79 C 7 y). The interval C > 1.9 corresponds to BR(H —>■ 77 ) > 3 
BR sm(H — > 77 ), while, excluding C Z 7 > 2.7, one excludes BR (H —> Z 7 ) > 9 BR sm(H —> Z 7 ). 
The corresponding BR bounds for the Higgs decay into a dark photon at the ILC with = 250 
GeV, and ~ 250 fb _1 of integrated luminosity, are about an order of magnitude looser than the 
FCC-ee ones. 

Of course, in order to fully assess the e + e _ collision potential, one would need an estimate 
of the corresponding LHC sensitivity to the H 77 and HZ 7 couplings. This would in particular 
allow to figure out whether, in case of negative findings at the forthcoming LHC runs, there 
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will indeed be some unexplored parameter space left that can be covered by searches in e + e“ 
collisions. Studies to evaluate the LHC sensitivity to the H —y 77 channel at Run 1 have 
just been started. A parton-level analysis leads to a bound on BR (H —y 77 ) of about half 
a percent at 95% C.L. from the 8 —TeV data set [TT] . corresponding to the exclusion of the 
range C > 1.6 (cf. Figure [d]). On the other hand, since most of the background to the 
H 77 signal in pp collisions comes from jet- and photon-mismeasurement effects, a thorough 
detector simulation would be needed to make a robust sensitivity statement on the present 
data set at the LHC. Studies at 13—14 TeV are expected to increase the sensitivity on a purely 
signal-statistics basis, but will be affected by harsher experimental conditions, that make the 
reconstruction of relatively small missing transverse energies quite critical. The extrapolation 
of present 8 —TeV results to larger pp c.m. energies will then require an even more careful 
analysis. Similar considerations apply to the LHC searches for a H —y Z 7 signal, where the 
expected missing transverse energy is even lower than in the H —> 77 case. 

Predictions for the BR (H — y 77) in the framework of the Flavor model proposed in [ 8 ] 
have also been presented. Due to non-decoupling effects, BR (H —» 77) turns out to be directly 
proportional to the mixing parameters in the scalar messenger sector, and to the U ( 1 )f 
coupling constant a, times some SM couplings. Remarkably, quite high £ q> i mixings are required 
in order to generate the Yukawa couplings radiatively, and avoid large fine tuning in the Higgs 
sector [Sj. At the same time, large (but still in the perturbative range) a couplings are expected, 
in order to avoid fine-tuning among the dark fermion charges. As a consequence, in the scenario 
|! 8 j, BR (H —> 77) can be naturally large, and in the ballpark of sensitivity of the LHC and 
future colliders. 

In conclusion, although realistic LHC sensitivity studies for the Higgs couplings to dark 
photons are still to come, one expects quite a lot of complementarity of the LHC and future 
e + e _ capabilities to probe new exotic Higgs-boson interactions with dark photons, thanks to 
the extremely clean e + e _ experimental conditions. 
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Appendix 


One of the most interesting property of the Flavor model in [ 8 ] is that the value of the U(l)p 
coupling strength in the dark sector, a, can be connected to the charge splitting among two 
dark-fermion generations in the lepton or quark sector. In the following, we present such a 
prediction. 

A crucial assumption in the model is the one of minimal flavor violation, implying that 
the only source of flavor violation comes from the U(l)p charges. We first define a (A) by 
normalizing to 1 the largest U(l)p charge, which is the one associated to the dark-fermion 
partner of the top quark. Then, according to Eq.fll7|), we get 


l/o;(A) ~ 


3 *4 

2vr 1 - q* Qi 



(33) 


where rn t is the top-quark mass, rn l stands for another generic quark mass, and q Q . is the 
U(l)p charge of the corresponding dark-fermion partner Then, by fixing the ratio of two 

q Qj charges in Eq. (33), one can predict d(A), as well as all the remaining U(l)p charges. For 


instance, by requiring a charge splitting of 10% among the U(l)p charges of the dark-fermion 
partners of the third quark generation, namely q U3 ~ 1 and q Dz ~ 0.9, we get b(A) ~ 0.13, and 
also q,; 2 ~ 0.87, q D2 ~ 0.82, and q Vl ~ 0.78, q Dl ~ 0.76, where for the different quark masses 
we assumed the central values in [50] . 

On the other hand, in order to obtain a weakly coupled U(1)f theory, with a (A) -C 10~ 2 , a 
charge splitting -C 1% is required, leading to an unnatural fine tuning among the U (1 )p charges. 
Remarkably, the same conclusion holds when extracting a (A) from the Eq. ( p3| ) applied to the 
purely EW dark sector. In this case, assuming q Ez = 1 and q E2 = 0.9 for the dark fermions 
associated to the r and /i leptons, respectively, one obtains a = 0.17, which is of the same order 
of the coupling strength for dark fermions in the quark sector. In order to avoid an unnatural 


fine-tuning among the U(l)p charges, Eq.(33) suggests a large, but maybe still perturbative, 
U(1)f coupling in the dark sector. 

Actually, in [ 8 ] the scale A can be many orders of magnitude above the messenger mass scale, 
and it is useful to have the a value at low energy (for instance, at the average messenger mass 
fh), as a function of the charge splitting of a pair of different dark fermions. Indeed, it is the 
low-energy a that enters BR (H —» 77 ) and BR (H —> 77 ). We then need first to connect a(A) 
to a(fh) by solving the appropriate renormalization group equations for the U(1)f /3-function. 
Due to the large energy gap between A and fh, all dark-fermion and messenger masses can be 
approximated to a common low-energy scale around fh, neglecting the running between different 
mass thresholds at low energy. Then, including the dark-fermion and messenger contributions 


to the one-loop beta-function, as well as the dark-fermion mass definition in Eq.(14), one can 
remove the explicit dependence on the A scale, by the expression 


a (A) = a(m) 


i + m 


4 + 3Ar c 2 
9 N c 


+ R l o- 
2 9 


(34) 


5 We have used the property that the loop function Y 0 (xi), as defined in Eq.( 18 ), has a weak dependence on 
the dark-fermion mass Mq. , and can be approximated to a constant. 
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Figure 12: Predictions for a at the low-energy messenger scale m, as a function of the U(l)p 
charge splitting in the third generation, defined by Qq Ds = q.Q v — S3). 


where R^' 1 are dehned in Eq.(20). The peculiar solution in Eq.(34) arises from reabsorbing the 
usual log(A/M 0i ) term into the dark-fermion mass dehnition in Eq.(14). 

In Figure [l2| the a behavior is shown as a function of the U (1 )p charge splitting £3 associated 
to the third generation of quarks, with qQ D3 = Qq U:i (1 — <^ 3 ) - We assume the charge normalization 
Qu 3 = 1 , and also qp 3 = 1 for the U(l)p charges of the dark fermions and messengers of the 
third-generation leptonic sector. The end point at £3 = 0.39 corresponds to a (A) = 1. 

By requiring a natural charge splitting that is not smaller than 20%, a turns out to be quite 
strong, 0.04 < a < 0.18, but still in the perturbative regime. Then, one obtains quite naturally 
large values for the H —» 77 branching ratio. 
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